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ABSTRACT 

Nonlinear  Langmuir  waves  in  a plasma  governed  by  tlte  dimension- 
less equations:  i9E/3t  = -V^E  + nE,  3^n/9t*  •=  V^(n  + g(|E|^))  are 
studied,  where  E is  the  complex  amplitude  of  the  high-frequency  elec- 
tric field;  n is  the  low  frequency  perturbation  in  the  ion  density 
from  its  constant  equilibrium  value;  and  g is  a given  function  of  |e(^. 
General  conditions  for  the  existence  or  nonexistence  of  a class  of 
multidimensional  solitary-wave  and  nonlinear  periodic  travelling-wave 
solutions  in  the  form  E(t,x)  = h(k*x-vt)  and  n(t,x)  = s(k«x-vt)  are 
established.  The  results  are  applied  to  the  special  cases:  (i)g(|E(^) 
“IeI^  corresponding  to  the  usual  pondermotlve  force,  and  (li)  g(|E|^)= 
K[l-exp(-|E|*)] , K is  a positive  constant,  representing  ion  density 
saturation. 
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1.  INTRODUCTION 

The  formation, interaction  and  collapse  of  nonlinear  Langmuir  waves  In 

1-10- 


plasmas  have  been  studied  extensively  In  recent  years?  In  most  of  the 
existing  works,  attention  is  focused  primarily  on  the  formation  and  inter- 
action of  solitary  waves.  Exact  expressions  for  these  solitary  waves  for 


various  regimes  have  been  obtained  only  for  the  one-dlroenslonal  case.  Re- 


cently^ Gibbons  et  al"'  discussed" the  possibility  of  existence  of  solitary 

, * < . • i . It/-.''' 

Langmuir  waves  for  higher  dimensions.  In  this  study, obtain  conditions 


10' 


for  the  existence  or  nonexistence  of  multidimensional .nonlinear  Langmuir 


travelling  waves,  including  the  periodic  and  the  usual  solitary  waves.  ^ 

1 


We  begin  with  the  following  basic  equations  describing  the  nonlinear 
Interaction  of  high-frequency  electron  oscillations  with  an  ion  fluid: 


i8E/9t  - -V^E  + nE, 


a^n/at^  - V^(n  + g(!E|*)), 


(1) 


where  1 ■/-I;  E » (Ej^, . . . ,Ej^)  is  the  complex  amplitude  of  the  high-frequency 
electric  field  6 given  by 


6(t  ,x)  " Re{E(t,x)  exp  (-lojpt)  } ; 


(2) 


n Is  a real  quantity  corresponding  to  the  low-frequency  perturbation  in  the 


ion  density  from  its  constant  equilibrium  value  n^;  and  g is  a specified 


real-valued  function  of  |e|^.  Here,  we  have  used  dimensionless  quantities. 


The 


units  of  time  t,  spatial  coordinates  x • (x^ , . . . ,Xj^) , electric  field 


and  density  n are,  respectively,  3q/(2o«i)p),  (3rp/2)(qa)  , SCqC'n^itT/ 3) 


and  Aqan  /3,  where  o is  the  electron-ion  mass  ratio  m /ro. ; q ■ T/T  ; T ■ 
o e 1 e 


T^  + T^;  T^,  Tj  the  electron  and  ion  temperatures  respectively;  r^  the 


electron  Debye  radius,  and  a)^  the  plasma  frequency.  The  function  g is 


r .3 

: 1 


1 


Introduced  here  so  as  Co  permit  Che  consideration  of  a wide  class  of  non- 
linear effects  such  as  satviratlon. 

2.  TRAVELLING  WAVE  SOLUTIONS 
N N 

Let  IR  and  C denote  Che  N-dimenslonal  real  and  complex  Euclidean 

spaces  respectively,  and  C^(  IR;V)  the  space  of  all  m-tlmes  continuously 

differentiable  functions  defined  on  R and  caking  their  values  in  the  vector 

space  V.  The  norms  for  and  are  denoted  by  ||  • ||  and  1*  1 respectively. 

N N 

The  dot  notation  is  used  to  denote  the  usual  scalar  product  on  R or  C . 

Let  k be  a specified  unit  vector  in  R^  and  v a given  real  number 
corresponding  to  a constant  dimensionless  velocity.  We  seek  travelling- 
wave  solutions  of  (1)  in  the  form: 


E(t,x)  " h(k*x  - vt) , 
n(t,x)  ■ s(k*x  - vt) , 


(3) 


N 

where  h and  s are  undetermined  functions  in  ^ re- 

spectively. For  physical  reasons,  we  shall  restrict  h and  s to  functions 
such  that  |h(^)|  and  |s(^)|  are  uniformly  bounded  on  IR,  where  ^ “ k*x  - vt. 
In  particular,  we  shall  consider  multidimensional  solitary-wave  solutions 
which  are  analogous  to  those  in  the  one-dimensional  case.  Here,  we  require 
that  |h(^;)|  and  tend  to  finite  values  as  |5|  -► 

Substituting  (3)  into  (1)  leads  directly  to  the  following  equations 


. . . ,h^)  and  s: 

-Iv  dh/d^  + d^h/dC*  - s(C)h(^), 

(4) 

(v^-  l)d*s/d£;^-  d^g(  |h|  *)/dC*  , 

(5) 

where  we  have  adopted  the  rectangular  Cartesian  coordinate  system. 


Equation  (5)  can  be  Integrated  to  give 


L 


(v^-l)s(S)  - g(lh(^)|2)  + + C, 


(6) 


where  C and  C are  Integration  constants.  From  the  boundedness  requirement, 
we  set  C ■ 0.  Assuming  that  v^fl,  we  can  solve  for  s(5)  in  (6)  and  substi- 
tute it  into  (4)  to  give  a complex  diiferentlal  equation  for  h: 


dVdC^-  ivdh/d^  - (v^-l)-ng(  |h(0  1*)  + C}h. 


(7) 


It  is  advantageous  to  rewrite  (7)  in  polar  form.  Let  hj(^) 


Aj(C) 


exp(i9j(^)),  j - 1,...,N.  Then,  we  have 


d^A^/dS"  + Aj0’(C)(v-0'a))  - (vM)-»Aj{g(||A|2)  + c}. 


d^0yd^^  - (v-20^(U)d(ln  Aj)/d£;,  j - 1,...,N, 


(8) 

(9) 


where  A - (Aj^, . . . ,Aj^)  , ||a||  - |h|  and  0j-d0j/d5. 
Equation  (9)  can  be  integrated  to  give 


e^(0  - (v-PjAj‘)/2. 


(10) 


where  - A^ (0) (v-20j (0))  . 


Substituting  (10)  into  (8)  leads  to  the  following  differential  equa- 


tions for  Aj : 


d^A^/d^^  - f(p^,C,A)Aj,  j - 1,....N. 


(11) 


where 


f(Uj,C,A)  - (vjAj"-  v2)/4  + (v^-l)-‘{g(||A|*)  + C}.  (12) 


Evidently,  given  0,0^ (0) ,Aj (0) ,Aj (0) , J «■  1,...,N,  (11)  can  be  inte- 


grated independently.  Since  f is  a function  of  C and  (depending  on 


\ t 


Aj(0)),  (11)  must  be  solved  with  Initial  conditions  at  ^-0  which  are  con- 
sistent with  the  Aj{0)  in  Vij  • Also,  only  those  portions  of  solutions  of 
(11)  with  A(^)i0  (i.e.  Aj(C)>0,  J-1,...,N)  are  meaningful  here. 

We  note  that  (11)  can  be  rewritten  in  the  form: 


d^Aj/d^z  . 3U/3Aj.  j-l....,N, 


(13) 


where 


U(A,p,C)  - U^(|lAJ^C)  - U^/(8A^). 

MV 

2L'j^(|a!|^,C)  -1  {<g(n)  +Y}dri, 


(14) 


(15) 


< - (v^-l)"‘  , y - (v^-D'^C  - vV4 


(16) 


and  y - (Pj^, . . . , Uj,) . A first  Integral  of  (13)  is  given  by 

N 

I(A(0.A'(0)  = llA’(f;)r-  2Uj(||A(C)f  ,C)  - ^ Pj^Aj2(C)/4  = C^,  (17) 


3-1 


where  |1a' (f;)|| 


(O/dQ^  and 


Cj  - ||a’(0)|12-  2Uj^(||a(0)P,C)  “ ^ (v  - 20j(O)). 

3-1 


(18) 


Evidently,  if  for  some  j,  then  I(A(5) ,A’ (5))— as  ||a'(C)|1  and 

||a(C)11-^0.  Since  is  finite  for  finite  |1  A' (0)||  ,1|a(0)|1  ,C  and  6 j (0) , 3= 
1,...,N,  therefore  there  do  not  exist  solutions  of  (13)  or  solitary-wave 
solutions  of  (7)  such  that  ||a(C)|1  and  ||A'(5)1|~*0  as  l^j— " when  Pjl^O  for 
some  J . 

In  what  follows,  we  shall  focus  attention  on  the  particular  case 
where  p * 0.  Here,  we  have 


4 


(19) 


(0  = 0^(0)  + vC/2, 

as  a solution  of  (10)  or  (9).  Note  that  =■  0 when  (0)»0  and/or  (0)“ 
v/2.  This  implies  that  along  any  trajectory  of  (8),  (9)  starting  from  a 
point  z(0)  * (A.(0)  ,A' (0) , 6(0) , 0’ (0) ) in  the  set  2 »{(A,A' ,0,6' ) 

A^(v-20p  = 0,  j=l,...,N},  its  corresponding  phase  0(5)  = (0^(5)  , . . . ,0jj(5)) 
has  the  form  (19).  In  this  case,  f no  longer  depends  on  A(0)  and  0’(O), 
and  (13)  reduces  to 

d^A^/dS^  = 3Uj^/3Aj,  j-l,...,N.  (20) 

2N 

The  equilibrium  points  of  (20)  are  points  (A^,0)  in  IR  such  that  A^  are 
the  stationary  points  of  or  the  roots  of  the  equation  f(0,C,A)A  * 0. 
Obviously,  A^s  include  A = 0 and  all  those  A's  satisfying  g(|lA J*)=v*(v^-1)/A 

— c . 

To  obtain  some  qualitative  Information  on  the  solutions  of  (20),  we 
first  derive  a differential  equation  for  u(5)  “ I1A(5)B^-  By  direct  com- 
putation, 

d^u/dP^  » 2|!a' (5)P+2A(5)-d2A/d5^  = 2|A'(5)P  + 2uf(C,u),  (21) 

where  f(C,||AP)  » f(0,C,A)  as  defined  by  (12).  Along  an  Integral  curve 
(17)  corresponding  to  a fixed  and  y*0,  (21)  can  be  rewritten  as 

d^u/dC^  . 2{uf(C,u)  + + 2Uj^(u,C)}  - P(u,C,Cj^).  (22) 

Its  solution  starting  with  Initial  conditions 

u(0)  - !lA(0)p,  u'(0)  - 2A(0)-A'(0)  (23) 


satisfying 


11a' (0)11 ' - + 2lI^(llA(0)||*.C)  >0 


(24) 


describes  the  evolution  of  ||a(Q  ||  with  along  the  Integral  curve. 
A first  Integral  of  (22)  Is  given  by 

u(C) 


(u'(0)^  - (u'(0))^  + 


L 


2P(n,c,c^)dn  ^ Q(u,c,Cj^,u’(o)), 


(25) 


where  u'  ■ du/dC.  Equation  (25)  is  valid  only  when  Its  right -hand-side 
Is  nonnegative.  An  Implicit  expression  for  ||A((i)||^  can  be  obtained  by 
Integrating  (25); 


jA(OP 

Q(n,c,Cj^,u'(o))  * dn  - . (26) 

l|A(0)P 


Note  that  If  an  explicit  expression  for  ||a(C)||^  Is  obtainable  from  (26), 
then  A(^)  can  be  determined  by  integrating  each  equation  In  (11)  Inde- 
pendently with  vi*0. 

In  the  sequel,  we  shall  establish  conditions  for  the  existence  or  non- 
existence of  solutions  of  (20)  having  the  property  that  ||a(£;)||— as  1^1— 
or  solitary-wave  solutions  of  (7)  with  U“0. 

Theorem  1:  If 


•cg(u)  + Y > 0 (27) 

for  all  u>0,  then  there  do  not  exist  solutions  of  (20)  such  that  )|A(0))|>0 
and  |(a(C)||— ►O  as  U|-*«. 

Proof;  Condition  (27)  is  equivalent  to  f*(C,u)>0  for  all  u>0.  In  view 
of  (21),  we  have  d^u/dC^X)  Implying  that  |1a(0|1*  along  any  solution  of  (20) 
Is  a convex  function  of  i.  Hence,  It  1?  impossible  to  have  ||A(0)  ||  > 0 and 


6 


I 


|a(C)  ||-K)  as  Ul-^- 1 

Note  that  for  the  subsonic  (v^<l)  and  supersonic  (v^>l}  cases,  (27) 
Implies  that  g(u)  la  uniformly  bounded  above  and  below  by  v*(v*-l)/4-C 
respectively.  Also,  If  (27)  Is  a strict  Inequality,  then  (A,A')"(0,0)  is 
the  only  equilibrium  point  of  (20). 

Theorem  2:  Assume  that  the  following  conditions  are  satisfied: 

(I)  v^(v^-l)  > 4C  and  v^<l; 

(II)  g is  a strictly  monotone  increasing  function  in  Cj^(lR,]R)  with 
g(0)«0  and  there  exists  a positive  number  Uj^<<»  such  that 

J.u 

I ^g(n)dn  “ {v* (v^-l)/4-C}u  (28) 

0 ^ 

and 

Jg(ri)dn  > {v^(v^-l)/4-C}u  forallu>u,  . (29) 

0 ^ 

Then  (20)  has  a solution  A(O>0  for  all  56®.  with  ||A(0)|1>0  and  ||A'(0)||"0 


such  that 


|1a(0|1  and  llA’(^)ll-K)  as  \^\ 


Proof : First,  we  note  from  (15)  and  (17)  with  ti>*0  that  for  a solution  to 
have  property  (30),  the  Initial  condition  (A(0),A'(0))  must  satisfy 


Cj  - |1a’(0)|12-  2U^(||A(0)p,C)  - 0. 


We  shall  show  that  under  condition  (1),  C^^-O  Implies  property  (30).  From 

(17),  it  is  evident  that  when  (A(0),A'(0))  satisfies  (31), its  corresponding 
trajectory  is  a zero-level  curve  of  I(A,A')  defined  by 


KA.A')  - ||A'|2-  2Uj(||a||*,C)  - 0, 


or  the  points  along  the  trajectory  belong  to  the  set 


7 


I-‘(0)  - {(A.A')6»^”:||aT-  2U^(|aP,C)}.  (33) 

Obviously,  the  equilibrium  point  (A,A' )«(0,0)6  I ”*  (0) . Now  the  foregoing 
implications  can  be  established  by  verifying  that  (0,0)  Is  the  only  equili- 
brium point  in  I~^(0),  moreover.  It  la  a saddle  point. 

Let  (A^,0)  be  an  equilibrium  point  of  ( 2C0  with  |a^|>0.  Then,  A^ 
must  satisfy  f(0,C,A^)”0  or 


ikJ 


I 

I 

) 

Under  condition  (1),  ( 3^^^/ 3A*)  |^_q>0,  so  (A,A’)*(0,0)  Is  a sad- 

dle point.  Moreover,  it  is  the  limit  point  of  some  trajectory  lying  in 
l''*(0)  as  1^|— w.  Hence,  implies  property  (30). 

Next,  we  must  verify  that  there  exist  points  (A(0) ,A’ (0))»(A{0) ,0) 
with  I A(0)9  >0  such  that  Cj^-0.  From  (15)  and  (31),  this  corresponds  to 
finding  a ||A(0)||>0  such  that 

Ja(o)P 

2UJ^(lA(0)|(^C)  ^ I {i«(n)  +y}  dn  - 0,  (39) 

■'o 

which  Invlew  of  (29),  can  be  rewritten  as 

>(0)1 » 

Ia(0)P  - W(|a(0)P)  2 4{vS(vM)-4C)  I 0(0)40-  (*0) 

0 


Evidently,  under  condition  (11),  the  mapping  W has  a nonzero  fixed  point 

I1a(0)P«». 

We  have  established  that  there  exist  points  (A(0) ,0)  in  I“*(0)  with 
||a(0)||>0.  Now,  we  must  show  that  for  such  a point,  there  exists  a tra- 
jectory curve  lying  in  I“M0)  which  joins  (A(0),0)  and  (0,0).  This  is 

assured  when  I~^(0)  is  compact.  It  is  straightforward  to  show  that  I~^(0) 


Che  compact  interval  . Thus,  the  boundedness  of  I (0)  is  established. 

Finally,  since  only  the  nonnegative  solutions  of  (20)  are  meaningful 
here,  it  remains  to  show  that  for  a point  (A(0) ,0)£ 1 (0)  with  A(0)>0  and 
Ha(0)||>0,  its  corresponding,  solution  is  nonnegative,  Chat  is,  A(^)>0 
for  all  This  is  Immediately  apparent  from  the  fact  that  I"V0)“ 

l^*(0)ul_''(0)  and  I"*(0)nI  -*(0)  - {(0,0)},  where  I’*  (0)- {(A.A* ) 6 I"’ (0)  : 

A»  and  1"'(0)  -{(A,A')6I  "*(0):A4)),  since  0^(0, C)-0.  | 

Remarks : (R-1)  Theorems  1 and  2 give  respectively  sufficient  condi- 
tions for  the  nonexistence  and  existence  of  multidimensional  solitary-wave 
solutions  of  (1)  which  are  directly  analogous  to  those  for  the  one-dimen- 
sional case.  From  (33),  it  is  evident  that  I~*(0)  is  symmetric  about  A"0 
and  A’“0.  Also,  I^'(O)  and  1”*(0)  are  symmetric  about  A'*0.  Thus .under 
the  conditions  of  Theorem  2,  the  trajectory  curves  in  the  (A,A' )-space  cor- 
responding to  the  solitary-wave  solutions  of  (7)  satisfying  (30)  have  si- 
milar properties,  and  they  have  the  form: 

hj(C)  • A^ (5)exp  {i(0j (0)+vC/2)}  , j“l,...,N.  (42) 

(R-2)  Along  any  solution  of  (20),  its  corresponding  density  s(5)  can 
be  found  directly  by  solving  (26)  for  ||A(^)j|^  or  |h(5)(^  and  substituting  the 

A 

result  into  (6)  with  C"0.  Complete  knowledge  of  the  solution  A(^)  is  not 
necessary  here. 

We  note  that  if  the  assumptions  in  Theorem  2 are  satisfied  and  there 
exists  a positive  number  r^  such  that 

icg(r')  + y - 0 or  g(rp  - v*(v*-l)/4  - C,  (43) 

then  (20)  ha.s  an  uncountably  infinite  number  of  nonisolated, nonzero  equili- 
brium points  (A^,0)  such  timt  A^  lies  on  the  sphere  {A^€  jA^|"r^}. 


Nov,  we  show  thst  there  exist  solutions  A(^)  of  (20)  in  some  neigh- 
borhood of  these  equilibrium  points  such  that  their  norms  are  periodic  func- 
tions of 

First,  we  rewrite  (22)  In  the  form 

d^u/df/  ■ 3V(u,Cj^)/3u,  (A4) 


where 


V(u,Cj)  - 2 


{n(icg(n)  y)  + + 2Uj^(p,c)}dn, 


(45) 


and  the  initial  conditions  u(0)*|A(0)|^  and  u' (0)“2A(0) • A' (0)  are  chosen 
such  that  condition  (24)  is  satisfied.  It  can  be  readily  verified  that  if 
we  set  given  by 


C 


0 

1 


{Kg(n) 


(46) 


then  Is  a stationary  point  of  V(*,C°),  or  (u,u')"(r|,0)  is  an  equili- 
brium point  of  (44).  For  condition  (24)  becomes 

S.A(0)|^ 

|A*(0)1|2  - I {Kg(n)  + y}ir)  > 0.  (47) 

■'r' 

e 

Under  condition  (1)  and  (ii)  of  Theorem  2,  we  have  icg(n)"*Y>0  for  0§ri<rg»and 

<g(ri)'*Y^0  sH  n>r^.  Evidently,  (47)  Is  satisfied  if  and  only  if  |A(0)|j 

■r  . Hence,  the  only  solution  to  (44)  with  Initial  condition  (u(0),u'(0)) 
e 

satisfying  (47)  Is  the  equilibrium  solution  (u(^) ,u' (5))“(r* ,0)  for  all 

Now,  we  consider  the  solutions  of  (44)  with  and  Initial  con- 

dition (u(0)  ,u' (0))-(|  A(0)(|  * ,2A(0)-A' (0))  satisfying  (24),  where  6Cj  Is  a 
small  perturbation  of  about  C°.  Let  'ig(Cj^)  denote  a stationary  point 
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ot  V(-,C^)  or  a root  of  the  equation 

u{Kg(u)  + yH  + 2Uj^(u,C)  ■ 0.  (48) 

For  we  can  write 

u (C,)  - r*  + 6u  . (49) 

e l e e 

Clearly,  under  the  assumptions  of  Theorem  2,  6u^  depends  continuously  on  5C^ 
and  |6u^l— *0  as  |6Cj^|— ♦O.  Also,  since 

3^V(u,Cj^)/3u^  - 4{icg(u)  + y}+  2Kug'(u),  (50) 

we  have  (3*V(u,C- )/3u^) I ■ 2<r  g'(r  )<0,  or  u -r*  Is  a relative  maximum 
i U"r  e e 6 e 

e 

point  of  V(*,C°).  In  fact,  since  2Kug' (u)<0  for  all  u>0  and  Kg(rg)+V*0 
there  exists  a positive  number  e such  that  for  each  6Cj^,  |5Cj^  (<c.  Its  cor- 
responding u^  (Cj^)«r^+6u^  is  a relative  maximum  point  of  V(*  ,Cj+6Cj^) . Con- 
sequently, for  any  fixed  1 6Cj^  |<e , (44)  has  periodic  solutions  In  some 
neighborhood  of  the  corresponding  equilibrium  point  (u,u' )"(r*+6u^,0) 

They  are  given  by  the  solution  of 


(u’(5))*/2  » V(u(£:),C,+6C.)  - V(u^,C°+8C,)  + (u’)V2  (51) 

i.  X oil  o 

with  u(0)“u  , where  the  Initial  point  (u(0) ,u' (0))*(u  ,u')  is  sufficiently 
o o o 

close  to  (u,u')»(r^+6u^,0).  In  particular  ^e  can  choose  u^-||A(0)||  ^>0  and 
u^*2A(0) -A' (0)  such  that  condition  (24)  given  by 

liA’(0)P-  C°  + 6Cj  + 2Uj^(|A(0)||JC) 

- 6C^  + I {Kg(n)  + Y>dn  ^ 0 (52) 

r * 
e 

Is  satisfied.  This  Is  possible  for  any  positive  5C^.  The  existence  of 
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solutions  of  (20)  in  some  neighborhood  of  Che  equilibrium  points  (A^,0) 
with  IjA^y-r^,  whose  norms  are  periodic  functions  of  ^ follows  from  the 
fact  that  |6Ug|— »0  as  |6Cj^l— *0.  The  foregoing  result  can  be  summarized 

as  a theorem. 

Theorem  3;  Assume  that  the  conditions  of  Theorem  2 are  satisfied,  and 
there  exists  a real  number  rg>0  satisfying  (43).  Then  there  exist  solu- 
tions A(^)  of  (20)  in  some  neighborhood  of  the  equilibrium  set  {(A,A')G  IR^^; 
||A|*r^,  A'*0}  such  that  their  norms  ||A(C)|  are  periodic  functions  of 

Note  that  In  the  multidimensional  case,  the  periodicity  of  u(5)*|A(5) J * 
generally  does  not  imply  the  periodicity  of  A(^).  Since  the  energy  density 
of  the  electric  field  is  proportional  to  (h(5)|^»  solutions  with  periodic 
lh(5) I represent  oscillatory  energy  densities.  Evidently,  from  (6)  (with 
C»0) , the  periodicity  of  s(C)  is  implied  by  that  of  |h(5)|.  Now,  we  give 
a simple  sufficient  condition  for  the  nonexistence  of  periodic  travelling 
waves  in  the  sense  that  |h(C)  | and  s(C)  are  periodic  in  5. 

Theorem  4;  Suppose  that  the  following  conditions  are  satisfied: 

(1)  g is  a real -valued  continuous  monotone  increasing  function  defined 
on  IR  such  that  g(0)-0; 

(il)  v^(v='-l)  < 4C  and  v*>l; 

(ill)  the  initial  conditions  A(0)  and  A*  (0)  satisfy  ||a(0)||>0  and  C^^O, 
where  is  defined  in  (31). 

Then,  the  norm  of  the  corresponding  solution  A(^)  of  (20)  is  nonperiodic 

in  C. 

Proof : Consider  (22)  given  explicitly  by: 

-U 

d*u/dC*  - 4yu  + 2Cj  + 2tc  |ug(u)  + j g(n)  dn  | , (53) 
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where  Y and  < are  as  In  (16).  Under  condition  (i),the  term  in  (53)  is 

nonnegative  for  u>0.  From  conditions  (ii)  and  (iii),  we  have  y>0  so  that 
d^u/d(|^>0  for  all  u>0.  Since  u(0)*||A(0)||  *>0,  u is  a nonzero  convex  function 
of  ^ which  cannot  be  periodic.  | 

Theorem  4 gives  a sufficient  condition  for  the  nonexistence  of  super- 
sonic periodic  travelling  waves.  In  the  subsonic  case  (v^<l),  the  condition 
v^(v^-l)<4C  implies  that  y*^0»  Thus,  under  condition  (i)  of  Theorem  4,  we 

have  d u/d^  40  for  all  u>0  when  Cj^40,  which  implies  the  nonexistence  of  sub- 
sonic periodic  travelling  waves.  However,  corresponds  to 

JA(0)P 

||a'(0)P  4 y1|a(0)P  -hcJ  g(n)dn,  (54) 

■'o 

whose  right-hand-side  is  nonpositive.  Thus,  this  condition  can  be  satisfied 
only  in  the  trivial  case  when  A(0)-0  and  A'(0)“0. 

3.  SPECIAL  CASES 

Now,  we  apply  the  results  in  Section  2 to  equation  (1)  with  particular 
forms  of  g arising  in  physical  situations. 

3.1  g(|Ej^)  » |e|*;  This  corresponds  to  the  case  with  the  usual  pondermo- 

tive  force.  Here,  as  given  by  (15)  has  the  explicit  form: 

2U^(i|AP,C)  - yIIaP  + k|a|''/2,  (55) 

where  y and  tc  are  as  in  (16).  A first  integral  of  (13)  is  given  by 

N 

|a'(C)P-  y|A(C)|*-  k|a(C)|V2  - C^.  (56) 

J-1 

When  A(0)  and  9(0)  are  chosen  such  that  y*0,  the  equation  for  Aj(^) 
given  by  (20)  reduces  to 
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(57) 


d^Aj/d^^  - (Y  + k|a|*)Aj,  j-1 N. 

and  the  equation  for  u(^)"|a(C) P given  by  (22)  becomes 
d^u/dC^  - 3ku*  + 4yu  + 2Cy 
A first  Integral  of  (53)  is  given  by 


(58) 


(u’(5))V2  - Ku’(^)  + 2yu*(0  + 2C^u(C)  + C^,  (59) 

where  C2  is  an  integration  constant.  By  restricting  the  right-hand-side  of 
(59)  to  be  nonnegative,  we  have  the  following  implicit  expression  for  u(C)J 


+ 2Yn^  + 2Cj^n  + C^}  ^dn  - /2C,  cem.  (60) 

(0) 

Now,  we  apply  Theorem  1 to  this  special  case.  Clearly,  for  v*<l, 
condition  (27)  cannot  be  satisfied.  But  for  v^>l,  (27)  is  satisfied  when 

v* (v*-l)/4"C.  Under  this  condition,  there  do  not  exist  solutions  of  (20) 
or  solitary-wave  solutions  of  (7)  such  that  |h(0)|>0  and  |h(C)|-K)  as  UI-. 
To  apply  Theorem  2 to  this  special  case,  we  see  that  under  condition  (1), 
g(u)-u  satisfies  condition  (11).  Thus,  under  condition  (i),  (57)  has  so- 
litary-wave solutions  satisfying  (30).  For  such  solutions,  we  see  from 
(56)  with  y»0  that  must  be  zero.  Also,  from  (59),  Cj  must  be  zero, 
since  u(0)-|| A(0)jl  ^>0  and  u*  (0)-0,| A(0)  | must  satisfy 


I 


|iA(0)p  - -2Y/tc, 

where  Y^O  and  k<0  under  condition  (i).  Thus,  (60)  becomes 

^iiA(or 

I {ten*  +2Yn^}’”dn  - /2C,  ^6®. 

^2y/1k1 


(61) 


(62) 
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It  follows  that 


l|A(r,)  H - (2>/ |k  |)  ^ sech(/YC) . C^IR, 


and  in  view  of  (b),  we  have 


s(Q  * (v^-1)"' {(2y/ |k:|)  sech^(/YC)  + C}. 


Substituting  (64)  into  (11)  with  yj*0  leads  to  a set  of  uncoupled 
equations  for  A^s  given  by 

d^Aj/d^^  “ (y  K:|lAp)Aj  - y{l  - 2 sech  ^( /yO  . 

j-l....,N,  (65) 

which  can  be  integrated  Independently  to  obtain  Aj(^).  The  foregoing 
results  are  consistent  with  those  for  the  one-dimensional  solitary  waves^^ 
Turning  now  to  the  periodic  travelling-waves,  we  observe  that  under 
condition  (i)  of  Theorem  2,  (57)  has  an  uncountably  infinite  number  of 
nonisolated  equilibrium  points  (A^,0)  such  that  ||A^||  ^••r|*»-Y/K>0.  Also, 
(58)  can  be  written  in  the  form  of  (44)  with 


V(u,Cj^)  “ u(k:u^  + 2yu  + 2Cj^) . 


If  we  set  (2k)  . then  u^-y/K  is  a relative  maximum  point 

of  V(*,C°).  Now,  we  consider  the  solutions  of  (58)  with  initial  condi- 
tions u(0)“lj  A(0)|]  ^ and  u' (0)”2A(0)  ^A' (0)  satisfying  condition  (24)  given 
explicitly  by 

||A’(0)|12-  4 YviCO)  + ku2(0)/2  ^ 0.  (67) 

Let  U(Cj^)  denote  the  set  of  all  u(0)>0  satisfying  (67)  for  a fixed  Cj^. 

It  can  be  readily  verified  that  under  condition  (i)  of  Theorem  2,  we  have 
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(68) 


X JL 

ti(C^)-{u(0):[Y-(Y^+2|K|C^)M/|<ku(0)<(Y+(Y^+2|K|C^)2]/|K|} 

for  0^j^^Y*/(2<) 

U(C^)-{u(0):0$u(0)^[Y+(Y^+2|<|c^)^]/1<|}  for  C^^O,  (69) 

and  U(Cj^)  is  empty  for  C^<>^/(2k).  Note  that  for  Cj^"Y^/(2<),  U(Cj^)  con- 
tains only  the  point  u(0)“-y/<.  Thus,  from  Theorem  3,  if  condition  (1) 
of  Theorem  2 are  satisfied,  then  there  exist  solutions  A(^)  of  (20)  in 
some  neighborhood  of  the  equilibrium  set  {(A,A')C]R  :|A||  •>-y/k,A'“0} 

such  that  their  norms  |1a(^)  1|  are  periodic  in  These  solution  curves  cor- 

respond to  (56)  with  U“0  and  satisfying  0>C^>Y^/(2k)  . When  is  set 
to  zero,  we  have  solitary-wave  solutions  such  that  |jA(C)  ||  and  ||a' (C)  |l“^ 
as  |c"| — ^ as  given  by  (63).  In  this  case,  (u,u')“(0,0)  is  a saddle  point 
of  (58)  with  C^«0. 

Fig.l  shows  the  trajectories  of  (58)  with  Y”1  and  k»-2  in  the  (u,u')- 
plane  for  various  values  of  and  u(0)  satisfying  (67).  Note  that  for 
Cj^'=‘0,  (u,u' )>*(0,0)  is  not  an  equilibrium  point  of  (58).  In  fact,  these 
solutions  pass  through  the  origin  and  they  are  periodic  functions  of  E,. 

Fig. 2 shows  the  behavior  of  the  trajectories  in  the  ^-domain. 

Finally,  for  the  supersonic  case  v^>l,  we  have  from  Theorem  4 that 
if  v^(v^-l)<4C  and  A(0),A'(0)  satisfy  |)a(0)S>0  and  or 

||A'(0)|2  i.Y||A(0)P+  ic||A(0)||''/2  > 0,  (70) 

then  the  norm  of  the  corresponding  solution  A(^)  of  (57)  is  nonperiodic 
in 

3.2  g( )E I ^)»K(1 -exp(- 1 E I ^) ) : This  form  of  g,  with  K being  a positive  con- 

12 

stant,  has  been  proposed  by  Wilcox  and  Wilcox  to  represent  ion  density 
saturation.  For  this  g,  is  given  by 
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2U^(I(aP,C)  - (kK  +Y)||aP-  <K{1  - exp(-||A||^)},  (71) 

and  a first  Integral  of  (13)  is  given  by 

N 

||A’(0|P-(KK+Y)||A(Q||VKK{l-exp(-|A(0|P))  - 2uj^Aj^(0/4  - C^.  (72) 

j-1 

The  equations  for  the  Aj(^)'s  with  i^O  corresponding  to  (20)  have  the  form: 

d^Aj/d^*  = (y  + <Kll  - expt-|!A(Qj!^)]}Aj.  (73) 

The  evolution  of  u(f,)  - l|A(f,)|l^  with  f,  along  an  integral  curve  of  (73) 
specified  by  Is  governed  by 


d^u/dC^  - 2l2(Y+<K)u  + Cj  - kKI 1 -(l-u)exp(-u) ] } , 


which  has  a first  integral  of  the  form: 


(u’(C))V2  - 2(Y+<K)u"(^)  + 2u({;)[C,  + kK (exp (-u(£; ))-!)]  + C-.  (75) 


where  is  an  integration  constant.  By  restricting  the  right-hand-side 
of  (75)  to  be  nonnegative,  we  can  integrate  (75)  to  give  an  implicit  expres- 
sion for  u(£.). 

To  apply  Theorem  1 to  this  case,  consider  inequality  (27)  given  expli- 
citly by 


icK(l-exp(-u))+  Y >0 


for  all  u > 0. 


This  condition  is  satisfied  when 


v*>  1 and  Y ^ 0 


v*<  1 and  kK  + Y ^ 0* 


Thus,  under  (77)  or  (78),  there  do  not  exist  solitary-wave  solutions  such 
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Moreover,  from  Theorem  4,  when 


that  |h(0)l  > 0 and  |h(^)|— *0  as  |^|— «o  . 

Y >0  and  (A(0),A'(0))  satisfies 

I|a'(0)P  > (K:K+Y)||A(0)f  - icK{l-exp(-|A(0)|)  *) } , (79) 

then  the  norm  of  the  corresponding  solution  A(^)  of  (73)  is  nonperiodic  in 
in  K. 

Now,  consider  condition  (29)  in  Theorem  2 which  requires  the  existence 
of  a Uj^>0  such  that 

K(uj^  + exp(-Uj^)-l]  ■ (80) 

and  for  all  u ^ 

Klu  + exp(-u)~l]  > (-y/k:)u  (81) 

This  condition  is  satisfied  if 

y/K  < 0.  (82) 

Thus,  if  < 1 and  y ^ 0,  then  the  hypotheses  of  Theorem  2 are  satisfied. 
Hence  (73)  has  solitary-wave  solutions  satisfying  (30).  For  such  a solu- 
tion with  j|A'(0)||"0,  we  have  from  (72)  with  and  y*0  that  ||A(0)  ||  must 

satisfy 

||a(0)F  - KK(<K+Y)-‘'i  1-exp (-11  A(O)P)},  (83) 

which  always  has  a solution  j|A(0)||*>0  if  v^<l  and  y>0. 

Next,  wo  observe  that  if 

k:K/(y-++:K)  > 1 , (84) 

then  (73)  has  an  uncountably  infinite  number  of  nonisolated  equilibrium 
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» 


» 
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% “ 11^1^  “ lntK»C/(Y-K:K)l  > 0.  (85) 

Note  that  if  then  (84)  implies  (82)  and 

( Y + KK)  < 0.  (86) 

We  shall  verify  that  under  the  conditions  of  Theorem  2,  (74)  has  periodic 
solutions  in  some  neighborhood  of  the  point  (u,u ' )- (u^ ,0) . 

First,  we  rewrite  (74)  in  the  form  of  (44)  with  V given  by 

V(C^,u)  - 2(Y-HcK)u^  + 2(Cj^-k:K)u  + 2kKu  exp(-u).  (87) 

If  we  set  given  by 

C°  - 2(Y+'«)lnI(Y+KK)/(KK)]  + tcK  - (y+kK)  [l-ln(KK/ (y+tcK) ) ] , (38) 

then  u^  given  by  (85)  is  a stationary  point  of  V(*,C°),  or  (u,u' )«-(u^,0) 
is  an  equilibrium  point  of  (74).  At  this  point,  O^V(u,C°)/3u^)  | _ ” 

U"Ug 

2u^(Y+'^‘K).  Thus,  under  condition  (86),  is  a relative  maximum  point 
of  V(*,C°).  Now,  we  consider  the  solutions  of  (74)  for  various  values  of 
in  some  nelghboihood  of  C°,  with  initial  conditions  u(0)"||a(0)  ||  ^ and 
u' (0)-2A(0) • A' (0)  satisfying  condition  (24)  given  by 

|a'(0)||^  “ ‘^l  (rt-KK)u(O)  - k:K[  1-exp (-u(0))l  >0.  (89) 

As  in  Section  3.1,  let  (((C^^)  denote  the  set  of  all  u(0)>0  satisfying  (89) 
for  a fixed  or 

U(C^)  - {u(0)>0:  C^+(y-HcK)u(0)  > KK[l-exp(-u(0)) ] } . (90) 

It  can  be  readily  verified  that  if  v^<l  and  condition  (84)  Is  satisfied, 
then  U(Cj)  is  empty  for  all  < C*,  and 


: j 
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(91) 


U(cJ)  - {u*(0)}, 

where 

u*(0)  = ln[KK/(Y+KK)],  (92) 

C*  - (Y+KK)ln[KK/(Y+<K)],  (93) 

where  u*(0)  corresponds  to  the  point  of  tangency  between  the  line  yj^(u)= 
C*(kK)”'+  [1+y(<K)~']u  and  the  curve  y2(u)”l  - exp(-u).  Also, 


U(C^)  = {u(0) : u ^ 

u(0) 

< S} 

for  C*<Cj^<0, 

(94) 

U(Cj^)  - {u(0):  0 4 

u(0) 

4 fi} 

for  Cj^>0, 

(95) 

where  i5  and  u with  u<u  are  the  two  distinct  positive  roots  of  the  equation 

+ (Y-HcK)u  - KK[l-exp(-u)]  * 0.  (96) 

Thus,  from  Theorem  3,  if  v^<l  and  (84)  is  satisfied,  then  there  exist  solu- 
tions A(^)  of  (73)  in  some  neighborhood  of  the  equilibrium  set  { (A,A' 

|aP=  In [ kK/ (y+kK) ] , A'“0}  such  that  their  norms  ||a(^)  ||  are  periodic  in 
These  solution  curves  correspond  to  (72)  with  vi»0  and  satisfying  0>Cj^>C*. 
The  trajectories  of  (74)  with  Y“1»k*~2  and  K=1  for  various  values  of  and 
u(0)  satisfying  (89)  are  shown  in  Fig. 3.  Their  corresponding  trajectories 
in  the  ^-domain  are  shown  in  Fig. 4. 

4.  CONCLUSION 

We  have  shown  that  under  mild  conditions  on  the  nonlinearity  g,  (1) 
has  multidimensional  solitary-wave  and  periodic  travelling-wave  solutions 
(E(C),n(C))  in  the  sense  that  |e(C)|  and  n(^)  tend  to  finite  values  as  |^| 

-»*,  and  they  are  periodic  functions  of  ^ respectively.  Along  these  solu- 
tions, the  phase  of  E(^)  is  an  affine  function  of  Moreover,  u(^)*|E(^)|^ 
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Fig. 3:  Trajectories  of  (74)  with  yl  *^“-2  in  the  (u,u')-plane  for  u(0) 
satisfying  (89)  and  fixed  values  of  C (curves  1-5  correspond  to 
C -0.2707,  0.0,-0.1522,-0.2642,-0.3069  respectively);  curve  2 is 
tne  solitary-wave  solution. 
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satisfies  a scalar  second-order  ordinary  differential  equation  whose  solu- 
tions have  properties  similar  to  those  In  the  one-dlmenslonal  case.  Al- 
though in  this  study,  we  have  treated  only  the  case  with  electrostatic  waves 
(l.e.  VxE  > 0),  the  same  approach  may  be  used  to  obtain  results  for 
electromagnetic  waves. 
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